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ABSTRACT 

Two  major  tasks  performed  during  the  report  period  were  investigation  of: 
(i)  crack  tip  ignition  phenomena  under  rapid  pressurization  and  (ii)  crack 
propagation  in  burning  solid  propellants. 

Ignition  of  AP-based  composite  solid  propellants  located  at  the  ti.p  of 
an  inert  crack  was  investigated  botn  experimentally  and  theoretically^  The 
ignition  process  was  observed  by  simultaneously  using  a  high-speed  000 

pictures/s)  camera  and  a  fast-response  photodiode  systam.  He  a-  flux  to  the 
propellant  surface  was  measured  with  a  thin-film  heat -flux  gage.  Effects  of 
pressurization  rate,  crack-gap  width,  and  igniter  flame  temperature  on  the  . 
ignition  process  were  studied  experimentally.  Results  indicate  that  the 
ignition -delay  time  decreases  and  the  heat  flux  to  the  propellant  surface 
increases  as  the  pressurization  rate  is  increased,  The  decrease  in  ignition 
delay  with  increasing  pressurization  rate  is  caused  by  enhanced  heat  feedback 
to  the  propellant  surface  at  higher  pressurization  rates.  This  augmentation 
in  heat  feedback  to  the  propellant  at  higher  pressurization  is  believed  to  be 
a  result  of  a  combination  of  the  following  mechanisms:  heating  due  to 
compression-wave  reflection  at  the  closed  end;  heat  release  due  to  burning  of 
unreacted  igniter  species  (or  particles)  near  the  tip,  behind  the  compression 
wa,fe;  and  enhanced  heat  transfer  due  to  recirculating  hot  gases  near  the  tip. 

In  the  theoretical  investigation  of  the  tip  ignition  event  a  comprehen¬ 
sive  model  for  ignition  of  AP-based  composite  solid  propellants  was  developed 
and  numerical  solutions  were  obtained^  The  analysis  simulates  the  ignition 
process  of  a  propellant  sample,  located  in  a  stagnation  region,  under  rapid 
pressure  loading  conditions.  Specific > features  considered  in  the  model  in¬ 
clude:  a)  detailed  chemical  kinetics  information  for  the  ignition  of  AP- 
fcased  composite  propellants;  b)  two-dimensional  (axisymmetric)  geometry  for 
the  composite  propellant;  and  c)  -apid  pressurization  of  the  gas  phase.  An 
implicit  finite  difference  scheme  was  used  to  solve  the  set  of  transient, 
second-order,  coupled,  inhomogeneous,  nonlinear,  governing  partial  differen¬ 
tial  equations.  Numerical  solutions  revealed  a  number  of  important  events 
occurring  during  the  ignition  sequence.  These  include:  igniter  gas  penetra¬ 
tion  to  the  region  near  the  saiqple  surface;  combustion  of  unbumed  species 
upon  arrival  of  compression  waves;  heat  transfer  to  the  propellant;  pyrolysis 
of  oxidizer  and  fuel;  and  gas-phase  reactions  leading  to  ignition.  The  model 
predicts  the  experimental  observation  that  the  ignition  delay  time  decreases 
as  the  pressurization  rate  is  increased.  Various  ignition  criteria  consid¬ 
ered  showed  the  same  trend  as  ejqjetimen tally  measured. 

Crack  propagation  in  a  burning  composite  solid  propellant  subjected  to 
rapid  pressurization  in  the  order  of  10  GPa/s  was  investigated  experimentally. 
A  visual  record  of  the  event  was  made  using  a  high-speed  movie  <  amera.  The 
transient  pressurization  process:  was  recorded  by  high-frequency  pressure, 
transducers  in  the  combustion  chamber .  ^Tne  effect  of  pres surizat ion  rate 
on  both  crack  propagation  velocity  and  time  variation  of  crack  shape  was 
studied.  Experimental  results  indicated  that  the  crack  velocity  increa  ;es 
as  the  pressurization  rate  is  raised^  The  measured  dependence  of  crack 
propagation  velocity  on  pressurizatTcJil  rate  was  found  to  be  stronger  than 
that  based  only  upon  viscoelastic  material  property  variation  at  different 
loading  rates.  A  characteristic  difference  of  crack  geometry  variation  during 
propagation  was  observed  at  different  pressurization  rates.  At  higher 
pressurization  rates,  the  geometric  transformation  of  the  crack  tip  was 
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ABSTRACT 

(Continued) 

preceded  by  emergence  of  a  fan  region  immediately  above  the  crack  tip.  It  is 
believed  that  the  fan  region  is  created  by  micro -structure  damage  in  the  propel 
lant  at  the  crack  tip,  which  allows  the  bright  combustion  gases  to  spread  ahead 
of  the  tip  region  in  a  radial  fashion.  The  actual  geometric  change,  a  trans¬ 
formation  of  the  crack  shape  from  a  smooth  triangular  contour  into  a  square 
contour  with  a  jagged  leading  edge,  is  believed  to  be  the  growth  of  small  frac¬ 
tures  into  major  crack  branches. 
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I .  INTRODUCTION 


This  report  summarizes  progress  made  during  the  period  September  1,  1981  to  Aug 
gust  31,  1982,  under  the  project  entitled  "Ignition  of  Solid  Propellants  and  Propa¬ 
gation  of  Burning  Propellant  Cracks"  (Contract  No.  N00014-79-C-0762) . 

The  overall  objective  of  this  investigation  is  to  achieve  better  understanding 
of  the  ignition  mechanism  under  transient  pressure  loading  conditions  and  to  investi 
gate  propagation  of  burning  solid  propellant  cracks.  These  fundamental  studies  are 
expected  to  help  in  providing  insight  into  convective  burning  and  evolution  of 
deflagration  to-detonation  transition  (DDT)  processes  in  damaged  solid  propellant 
grains.  Specific  objectives  of  this  study  are: 

1.  To  numerically  solve  the  theoretical  model1  for  ignition  of  AP-based 
composite  solid  propellants  under  rapid  pressure  loading  conditions. 

2.  To  study  the  effect  of  such  parameters  as  pressurization  rate  and 
oxidizer  particle  size  on  the  ignition  piocess. 

3.  To  validate  the  theoretical  model  by  comparing  predicted  results 
with  experimental  data. 

4.  To  examine  the  effect  of  various  ignition  criteria  on  predicted 
results. 

5.  To  observe  crack  propagation  process  in  burning  solid-propellant 
samples  under  rapid  gas  loading  conditions. 

To  measure  crack  propagation  velocities  in  burning  propellant  samples. 

7.  To  study  the  effect  of  pressurization  rate  on  crack  propagation 
velc -ity. 
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by  A.  K.  Kulkami,  M.  Kumar,  and  K.  K.  Kuo 

IT.  INVESTIGATION  OF  COMPOSITE  PROPELLANT  IGNITION 
2.1  Brief  Description  of  Ignition  Model 

The  modified  mathematical  model  is  presented  briefly  below.  More  details  of 
the  model  can  be  found  in  our  last  annual  report.*  The  model  simulates  ignition 
of  a  composite  propellant  in  a  stagnation  region.  The  physical  model  considers  a 
cylindrical  oxidizer  particle  having  radius  R^  and  length  L  embedded  in  a  cylin¬ 
drical  binder  of  radius  R2  *-n  suc^  a  way  that  initially  the  oxidizer  and  ^'lel 
binder  surfaces  are  in  the  same  plane  (see  Figs..  1  and  2).  For  the  coordinates 
shown  in  Fig.  2,  the  solid-phase  equations  are 
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Fig.  1  Typical  Composite  Propellant  Structure 
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Fig.  2  Cross-Sectional  View  of  a  Statistically 
Averaged  Element  and  Finite  Difference 
Grid  Pattern 
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The  expression  for  qp' s  is  similar  to  that  of  q”^  ; 
The  gas-phase  conservation  equations  are 
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where  j  *  1,  2,  iand  4  fcr  the  gas  phase  species,  oxidizer,  NHj,  HClQ^and  fuel 
respectively. 

The  gas-phase  momentum  equations  are  replaced  by  the  measured 


P  «  P(t) 


The  equation  of  state  for  the  gas  phase  is 


PW 

P*  =  V 

The  following  basic  assumptions  are  employed  in  the  mathematical  model; 
1)  The  solid  propellant  and  gas  phase  are  t'.’o- dimensional  axisymmetric. 


s 


2)  The  gas-phase  pressure  is  a  prescribed  function  of  tine,  and  is  uniform 
in  a  small  region  adjacent  to  the  solid  propellant. 

3)  Gases  obey  the  perfect  gas  law. 

4)  Chemical  reactions  and  pyrolysis  processes  can  be  described  by 
Arrhenius  expressions. 

5)  Binary  diffusion  coefficients  of  all  species  are  equal. 

6)  Prior  to  ignition,  displacement  of  the  propellant  surface  due  to 
regression  can  be  neglected. 

7)  The  radial  bulk  velocity  is  much  smaller  than  the  axial. 

The  initial  and  boundary  conditions  for  the  solid  phase  (both  fuel  and  oxidizer) 
are 

at  t  -  0  :  T(0,r,t)  ■  T. 

(11) 


at  x  ■*>  -  •  j  T  *  T. 


r  -R  + 


*1 


_  and 


a  -  -L" 


-  -L+ 


on  r  -  0  s  “  0 
dr 


and  on  r  ■  R„ 


3T 

3?"  0 


at  *  -  * 


a-g 


k  s. 


3T 


•  M  "  k8  3* 


>  +  +  Vgc  + 

l-g  8 


rfc  P,T(<V  -  c  )  +  dM 

b8  8  Pg  P8  •- 


8 


at  z  -  -L 


3T 


.  t,  Oi 

, s  3z 


-  k  II 

-L  Ox, g  3s 


-L+  + 


‘WW5* 


-6  «• 
Ox 


(12) 

(13) 


(14) 


(15) 


(16) 


(17) 


6 


9 


at  r 


•*X 


id  -L  <  a  <  0 


3T 


•  -  k_  -Si  .  a 

S.a  3r|  5^+ 


_  .  3T 
“Oa.a  "Sr  1^“ 


(18) 


Initial  and  boundary  conditions  for  the  gas-phase  equations  (for  j  »  1,  2,  3,  4) 

are 
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The  following  reactions  are  considered  in  the  model. 
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The  reaction  rate  constant  is  approximated  by 


k  •  2  x*e”*^*u* 

If  the  gas-phase  mass  fractions,  Y.,  Y2,  Yj,  Y^,  Yj, 
ammonia,  perchloric  acid,  fuel,  and  products,  respectively 
be  expressed  as  follows. 
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The  following  paragraphs  list  differences  in  the  model  presented  above  fro a 
that  given  in  Ref.  1.  (Complete  solution  of  the  ignition  aodel  for  composite  pro¬ 
pellants  is  presented  in  Section  2.3.) 

1.  In  our  previous  version,  axial  diffusion  of  species  in  the  gas  phase  was 
neglected;  it  has  been  reinstated  in  the  current  aodel.  This  was  done  because  it 
was  found  that  very  near  the  propellant  surface  the  axial  species  concentration 
gradients  can  be  extreaely  steep  which  are  cougar able  to  the  axial  convection  or 
radial  diffusion. 

2.  In  the  sol id -phase  energy  equation,  the  inhomogeneous  terms  q!"  and 

vX  p  S 

q"'  included  heat  generation  due  to  photochemical  reactions.  To  model  this  term, 
r  ,s 

values  of  Qpc^  (rate  of  heat  generation  due  to  photochemical  process  per  unit  wave¬ 
length  of  radiation)  and  (efficiency  for  photochemical  process  at  wavelength  X) 
are  needed  for  both  the  fuel  and  oxidizer.  These  are  not  available  in  the  literature, 
and  also  the  contribution  by  the  photochemical  reaction  was  thought  to  be  small  in 
AP/PBAA  propellant.  Therefore,  this  term  is  neglected  in  the  current  model. 

3.  The  total  energy  loss  by  radiation  from  the  propellant  surface  to  the 
environment,  cs^s>  (which  appears  in  the  energy  balance  at  the  solid-gas  interface 
was  estimated  to  be  small  compared  to  convective  and  conductive  heat  transfer 
during  the  ignition  process,  therefore  it  has  been  neglected. 

4.  Heat  release  due  to  reactions  at  the  solid-gas  interface,  q"  ,  was 

5-g 

restricted  to  latent  heat  of  gasification  represented  by  r.  p.[Ah-  -  Ah*  ].  No 
heterogeneous  chemical  reactions  were  considered  due  to  the  lack  of  adequate  infor¬ 
mation  on  heterogeneous  chemical  reactions  in  AP/PBAA  propellant  ignition.  Heat 
generation  due  to  reactions  at  the  fuel-oxidizer  interface  in  the  solid  phase  is 
also  neglected  because  previous  investigations  did  not  find  evidence  of  such 
reactions  for  AP  composite  propellants. 


5.  Far  avry  from  the  propellant  surface  (z  •+■*),&  new  boundary  condition 

had  to  be  added  because  of  the  inclusion  of  axial  diffusion  term  in  the  species  equa- 
3Y 

tion.  Therefore,  »  0,  was  imposed  to  prevent  any  species  diffusion  at  the 

o  Z 

far-field  boundary.  It  should  be  noted,  however,  the  species  can  still  be  trans¬ 
ported  across  this  boundary  hy  convection. 

3T 

6.  Far  away  from  the  propellant  surface  (z  **“  ®)»  3  because  the  spatial 

variation  in  temperature  is  not  expected  to  be  present  there. 

7.  There  were  some  minor  typographical  errors  in  the  model  presented  in  the 

1  •  • 
last  annual  report:  a)  the  term  (1/T)  in  the  expressions  for  to^'and  co^'were 

replaced  by  T.  b)  The  negative  sign  appearing  in  the  equation  for  q"' should  be 

positive,  c)  Some  signs  in  the  energy  balance  at  the  interface  z  *  -  L  have  also 

been  corrected. 

2 . 2  Coaputt r  PiogTua  Development 

A  finite  difference  technique  has  been  used  to  develop  a  computer  program  to 
solve  the  composite  propellant  ignition  model.  A  set  of  nine  coupled,  nonlinear, 
inhomogeneous,  partial  differential  equations  along  with  a  complete  set  of  initial 
and  boundary  conditions  wake  up  the  mathematical  model.  As  a  issult  of  this 
cosplexity,  no  analytical  solutions  exist  or  could  be  developed;  therefore,  numerical 
solution  was  the  only  alternative  to  obtain  theoretical  predictions. 

2.2.1  Special  Features  of  the  Numerical  Scheme 

Several  common  features  described  in  the  following  were  used  to  obtain  an 
accurate  solution.  Because  large  temperature  and  species  concentration  gradients 
exist  near  the  propellant  surface,  a  coordinate  transformation  was  performed  on 
the  z-axis  to  obtain  an  exponentially  finer  grid  spacing  in  the  real  coordinate 
while  maintaining  a  constant  grid  spacing  in  the  transformed  coordinate. 

Solid  phase:  e  -  «xp(A#s)  -  1  or,  x  -  ^  in  (14s) 
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In  order  to  control  possible  stability  problems ,  a  generalized  Crank-Nicolson 
technique  was  employed  to  solve  the  system.  Allen's  method  was  utilized  to 
approximate  the  radial  derivatives  with  a  variable  grid  size  while  a  three  point 
central  difference  scheme  was  used  to  approximate  the  axial  derivatives  with  a 
constant  grid  spacing. 

Quasi linearization  of  the  inhomogeneous  terms  was  required  because  the  source 
terms  in  the  governing  equations  are  implicitly  dependent  upon  temperature  and 
species  concentration.  Finally,  a  successive  overrelaxation  iterative  scheme  was 
chosen  to  solve  a  larger  set  of  algebraic,  finite-difference  equations. 

Appendix  A  describes  the  details  of  the  abovementioned  numerical  techniques. 

2.2.2  Methods  Used  in  Checking  Numerical  Solutions 

There  were  three  major  checks  made  on  the  computer  solution  to  ensure  that 
numerical  solutions  obtained  are  the  genuine  solutions  of  the  theoretical  model: 

i)  Comparison  with  available  analytical  solution  -  to  show  that  the  numerical 
solutions  for  simplified  geometric  configuration  under  known  surface 
heat  flux  agree  with  the  classical  analytical  solutions, 

ii)  Global  energy  balance  -  to  sum  the  net  heat  addition  to  all  elemental 

volumes  in  the  solid  phase  and  to  verify  that  the  overall  energy  balance 
is  achieved. 

iii)  Local  energy  balance  check  -  to  examine  whether  the  numerical  solution 

can  satisfy  the  governing  differential  equations  and  boundary  conditions. 

All  three  independent  checks  verify  that  the  numerical  solution  obtained  from 
the  computer  code  represents  the  correct  solution  of  the  partial  differential  equa¬ 
tions  in  the  theoretical  model. 


2.2.3  Current  Status  of  the  Computer  Code 

One  year  ago,  the  computer  code  was  only  capable  of  solving  simplified  cases 
in  which  measured  heat  flux  to  the  propellant  surface  was  used  as  an  input  to  the 
solid-phase  energy  equations.  Only  the  solid-phase  portion  of  the  whole  physical 
model  was  implemented  and  fully  debugged;  the  gas-phase  portion  was  partly  coded. 
Currently,  the  complete  model  including  both  gas-  and  solid-phase  equations  has  been 
implemented,  debugged,  and  tested.  A  layout  of  the  flow  diagram  of  the  computer 
code  is  given  in  Fig.  3. 

At  the  present  time,  the  computer  program  is  capable  of  predicting  ignition 
process  for  AP/PBAA  propellants  only  because  the  current  reaction  mechanism  is 
specific  to  the  above  propellant.  However,  the  program  can  be  extended  to  solve 
for  any  other  propellant  as  long  as  the  reaction  mechanism  for  that  propellant  is 
prescribed.  In  the  near  future,  the  ignition  mechanism  for  nitramine  propellants 
will  be  incorporated  into  the  model  and  will  subsequently  be  implemented  in  the 
computer  code.  Also,  the  computer  code  which  is  presently  executing  on  IBM  370/ 

3018  will  be  installed  on  PDP  11/23  to  reduce  the  computer  cost. 

2.3  Numerical  Solutions  and  Comparison  with  Experimental  Results 

2.3.1  Input  Parameters  and  Initial  Conditions 

Numerical  values  of  the  important  input  parameters  to  the  computer  programs 
are  listed  in  Table  1,  Thermal  properties  of  AP  were  taken  from  Ref.  8.  Thermal 
properties  of  PBAA  are  approximate  and  axe  identical  to  those  used  by  Varney  and 
Strahle.7  Thermal  conductivity  of  PBAA  was  deduced  from  the  average  density  of 
fuel  and  the  value  of  thermal  diffusivity  used  in  Ref.  7  for  fuel  binders.  Thermal 
properties  for  the  gas  phase  were  obtained  by  averaging  the  thermal  properties  of 
the  combustion  product  species  of  the  AP/PBAA  propellant.  The  thermal  conductivity 
is  assumed  to  vary  as  T0,7^.  The  composition  of  these  product  species  was  deter¬ 
mined  from  the  NASA-Lewis  chemical  equilibrium  calculation  program.10  Molecular 


weights  of  the  gas-phase  species  were  obtained  in  a  similar  fashion. 
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Numerical  values  of  chemical  kinetic  constants  are  given  in  Table  2.  Some 
of  the  activation  energies  given  in  Table  2  are  approximate  due  to  lack  of  kinetic 
data.  It  should  be  noted  that  chemical  kinetics  of  ignition  may  differ  signifi¬ 
cantly  from  those  of  steady-state  combustion.  It  is  generally  believed  that  AP 

2 

pyrolysis  consists  of  70%  exothermic  degradation  and  30%  dissociative  sublimation. 

The  activation  energy  of  this  decomposition  is  estimated  to  be  between  20-30  kcal/ 
2-7 

mole.  In  the  present  study,  the  value  was  taken  as  24  kcal/mole.  The  activation 

3 

energy  of  the  premixed  reaction  between  NH3  and  HCIO^  is  15.5  kcal/mole.  The 

pyrolysis  of  fuel  (PBAA)  is  endothexmic.  The  activation  f nergy  for  fuel  pyrolysis 

was  taken  as  113  kJ/mole  (27  kcal/mole) ;  Ninan  and  Krishman11  have  shown  that  for 

many  polybutadiene  binders  the  activation  energy  for  thermal  decomposition  is 

g 

around  110  kJ/mole.  (It  may  be  noted  that  Varney  and  Strahle  reported  the  activa¬ 
tion  energy  for  PBAA  decomposition  as  34  kcal/mole) .  The  pyrolysis  products  of 
fuel  and  oxidizer  react  in  the  gas  phase  to  form  the  final  products,  however, 
no  reliable  estimate  of  kinetic  parameters  for  this  reaction  is  available. 

The  initial  conditions  for  computations  were  consistent  with  those  encountered 
in  the  experiments.  Initial  pressure  and  temperature  were  0.1  MPa  and  300  K, 
respectively.  In  order  to  simulate  the  bright  reaction  zone  observed  near  the 
tip  (see  Refs.  1  and  2),  a  mixture  of  gaseous  oxidizers  and  fuel  species  was  assumed 
to  be  present  at  t  =  0.  Pressure-time  traces,  used  as  input  to  the  computer  pro¬ 
gram,  were  either  experimentally  measured  or  prescribed  for  a  given  pressurization 
condition . 

2.3.2  Ignition  Criteria 

Since  there  is  no  universally  accepted  ignition  criterion  (IC) ,  several 


ignition  criteria  were  employed  to  test  the  predictive  capability  of  the  model. 
The  onset  of  ignition  is  defined  as  the  time  (1)  when  the  local  maximum  of 
volumetric  gas-phase  heat  release  reaches  a  critical  value  (q"’  >c) ,  (2)  when 

g  >  IT.  3- A 

the  AP  pyrolysis  rate  near  the  radial  oxidizer-fuel  interface  reaches  a  critical 
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value,  (3)  when  the  axial  gas  velocity  at  r  =  1^  and  a  distance  away  front  the 
surface  is  zero  (i.e.,  when  the  surface  blowing  is  sufficiently  large  to  overcome 
the  incoming  gas  velocity  induced  by  pressurization  of  the  cavity),  (4)  when 
the  rate  of  maximum  volumetric  gas -phase  heat  release  reaches  a  critical  value, 
or  (5)  when  the  rate  of  AP  pyrolysis  rate  reaches  a  critical  value. 

2.3.3  Discussion  of  Predicted  Results 

The  general  sequence  of  events  leading  to  ignition,  as  revealed  by  the  numerical 
solution  is  as  follows. 

As  pressurization  of  the  crack  cavity  begins,  temperature  in  the  crack  increases; 
and  gases  flow  toward  the  crack  tip  at  high  velocities  (see  Fig.  4a).  The  rise  in 
gas  temperature  produces  a  temperature  gradient  at  the  solid-gas  interface,  heat 
is  transferred  to  the  solid  propellant,  and  heating  of  the  solid  propellant  begins. 

At  the  same  time,  as  the  gas  temperature  rises,  the  initial  unreacted  oxidizer  and 
fuel  species  (from  the  igniter  system)  react,  further  increasing  the  local  gas 
temperature  as  well  as  the  heat  flux  to  the  prope1lant  (see  Fig.  5).  This  results 
in  a  continuous  increase  in  the  propellant  surface  temperature  during  the  initial 
period. 

As  time  progresses,  the  initial  unreacted  species  present  in  the  gas  phase 
are  completely  consumed,  and  the  temperature  of  the  propellant  surface  continues 
to  rir«.  This  causes  slight  attenuation  in  the  heat  flux  to  the  propellant  sur¬ 
face.  The  heat  flux  increases  very  rapidly  during  the  initial  period  and  then 
levels  off.  The  axial  temperature  distribution  at  the  oxidizer  fuel  interface 
(r  -  Rj)  at  various  times  is  shown  in  Fig.  6.  As  seen  in  Fig.  6,  the  thermal  wave 
penetration  into  the  solid  increases  with  time.  Temperatures  far  from  the  sur¬ 
face  are  quite  uniform  (not  shown  in  the  figure) .  The  reversal  of  maximum  temp¬ 
eratures  with  time  at  20  ym  into  the  gas  phase  is  caused  by  a  relaxation  from 
the  initial  reactions  and  a  decreasing  pressurization  rate  with  time. 


Fig.  4  Schematic  of  Sequential  Variation  of  Axial  Velocity 
Field  Adjacent  to  Propellant  Surface 
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Fig.  6  Axial  Temperature  Variation  at  AP/PBAA  Interface 

(r=Rj)  for  Various  Times  (Average  dP/dt  =  120  GPA/s) 
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Since  the  thermal  properties  of  AP  and  PBAA  are  very  different,  there  is  a 
steep  and  discontinuous  temperature  gradient  near  the  AP/PBAA  interface.  Figure 
7  shows  the  radial  variation  of  propellant  surface  temperature  as  a  function  of 
time.  Because  PBAA  has  a  lower  thermal  diffusivity  and  conductivity,  its  surface 
temperature  rises  much  faster  than  that  of  AP.  The  gas-phase  temperature  near 
the  propellant  surface  also  changes  with  radius  and  is  strongly  coupled  to  the  solid- 
gas  interface  temperature.  Farther  away  from  the  propellant  surface,  the  radial 
gas  temperature  variation  diminishes.  Contrary  to  the  numerical  results  obtained 
here,  which  shows  that  the  surface  temperature  of  PBAA  is  much  higher  than  the 
AP,  Price  et  al.*1  deduced  from  studying  AP/PBAA  sandwiches  in  steady-state  burn¬ 
ing  that  the  fuel  surface  temperature  is  lower  than  the  AP.  If  the  thermal  dif¬ 
fusivity  and  conductivity  of  PBAA  were  larger  than  that  of  AP,  Price's  trend 
for  surface  temperature  would  be  observed;  however,  the  functional  relationship 
of  the  predicted  ignition  delays  versus  pressurization  rate,  etc.,  will  remain 
unchanged . 

Figure  8  shows  time  variation  of  maximum  volumetric  gas-phase  heat  release. 

The  gas-phase  reaction  between  the  unbumed  oxidizer  and  fuel  species  lasts  for 
about  10  ys.  During  this  interval  the  heat  release  rate  is  extremely  high  (the 
values  are  off  the  scale  of  Fig.  8) .  High  initial  concentration  of  the  species 
causes  the  rate  of  heat  generation  to  increase  abruptly  as  the  local  temperature 
increases  due  to  pressurization.  The  heat  generation  rate  peaks  and  then  drops 
to  almost  zero  (even  though  the  temperature  continues  to  increase)  since  the 
reacting  species  are  consumed  rapidly.  As  time  progresses,  gas  velocities  decrease 
because  of  pressure  rise  in  the  cavity  (see  Fig.  4b). 

Following  a  relatively  long  period  of  inert  heating,  the  AP  and  PBAA  start  to 
pyrolyze  when  the  propellant  surface  temperature  becomes  sufficiently  high.  The 
gasified  species  of  oxidizer,  fuel,  NHj  and  HCIO^  diffuse  away  from  the  surface. 


Radial  Distance,  r,  ym 


Fig.  7  Radial  Variation  of  Surface  Temperature  at 
Various  Times  (Average  dP/dt  *  120  GPa/s) 
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both  axially  and  radially.  As  shown  in  Figs.  9  and  10,  the  concentration  of  fuel 
species  is  highest  above  the  surface  at  r  ■  R2>  where  fuel  surface  temperature  is 
highest  (see  Fig.  7) .  The  fuel  mass  fraction  decreases  away  from  that  location. 

The  oxidizer  concentration  is  highest  at  the  AP-PBAA  interface  (r  *  R^),  as 
shown  in  Figs.  9  and  10,  where  the  AP  surface  temperature  is  highest  (see  Fig.  7). 

As  time  progresses,  the  mixing  of  fuel  and  oxidizer  species  continues  because  of 
diffusion  and  convection.  Since  local  temperatures  are  not  high  enough  to  support 
vigorous  gas-phase  reactions,  the  concentration  of  these  species  continues  to 
increase  as  a  result  of  accumulation. 

As  the  pressurization  process  continues,  the  temperature  of  the  gas-phase 
continues  to  increase.  At  the  same  time,  there  is  an  accumulation  of  fuel  and 
oxidizer  species  due  to  pyrolysis  of  the  propellant.  When  local  concentrations 
and  temperature  reach  suitable  values,  vigorous  gas-phase  reaction  between  pyTolyzed 
oxidizer  and  fuel  species  ensues  (see  Fig.  8).  Ignition  is  defined  as  the  attain¬ 
ment  of  high  gas-phase  reaction  rates  (see  Fig.  8)  and  simultaneous  high  surface 
pyrolysis  rates  (see  Fig.  4c  and  4d). 

2.3.4  Parametric  Studies 

The  critical  values  used  in  ignition  criteria  (see  Sec.  2.3.2)  were  obtained 

by  matching  the  predicted  ignition  delay  with  experimental  data  for  one  case. 

Comparison  between  predicted  ignition  delay  (for  first  three  ignition  criteria 

mentioned  earlier)  is  shown  in  Fig.  11.  It  is  cleaz  that  all  three  criteria 

12 

correctly  predict  the  experimentally  observed  trend  that  ignition  delay  time  is 
lower  for  higher  pressurization  rate.  The  last  two  criteria  also  show  the  same 
trend. 

Predicted  ignition  delay  time  shows  a  greater  dependence  on  the  choice  of 
ignition  criteria  at  lower  pressurization  rates.  The  criterion  based  on  zero 
velocity  (IC  3)  or  critical  burning  rate  (IC  2)  results  in  a  lower  ignition  delay 
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Fig.  10  Axial  Distribution  of  Oxidizer  (Y,)  and 
Fuel  (Y.)  Species  above  AP/PBAA  Interface 
(r-Rj)  (Average  dP/dt  *  120  GPa/s) 
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time  than  that  based  on  gas-phase  heat  release  (IC  1).  This  is  reasonable  because 
at  a  lower  pressurization  rate,  the  gas-phase  temperature  increases  at  a  slower  rate; 
therefore,  the  time  lag  botween  significant  pyrolysis  and  the  attainment  of  the 
critical  gas-phase  reaction  rate  is  longer.  Theoretical  results  overestimate  the 
ignition  delay  time  at  lower  pressurization  rates  and  slightly  underestimates 
the  ignition  delay  time  at  higher  pressurization  rates  (see  Pig.  11);  this  may  be 
caused  by  approximate  values  of  chemical  kinetics  constants  used  in  the  computation. 
Overall,  predicted  results  are  in  reasonable  agreement  with  the  measured  data. 

Another  possible  ignition  criterion  that  can  be  used,  but  was  not  employed 
for  results  presented  in  this  report ,  is  the  appearance  of  an  abrupt  rise  (spike) 
of  the  surface  tenperature  (see  Pig.  12)  and  the  gas-phase  temperature  adjacent 
to  the  radial  oxidizer-fuel  interface.  The  spike  occurs  because  of  exothermic 
reactions  on  the  AP  side  of  the  interface  and  endothermic  reaction  on  the  PBAA 
side.  This  spike  results  in  both  high  pyrolysis  rate  at  r  «  as  well  as  pronounced 
gas-phase  reaction  rate  near  the  circumference  of  the  interface  region.  As  evident 
from  Fig.  12,  the  spike  occurs  at  a  later  time  than  the  onset  of  ignition  based 
upon  previously  mentioned  ignition  criteria.  However,  the  predicted  trend  of 
ignition  delay  will  not  be  influenced  by  the  selection  of  this  ignition  criterion. 

It  should  be  noted  the  predicted  surface  temperature  at  the  onset  of  ignition 
is  about  450  K  which  is  considerably  lower  than  that  observed  for  steady-state 
burning.  This  is  reasonable  because  of  the  time  lag  between  the  onset  of  ignition 
and  steady  state  burning.  If  the  program  is  allowed  to  run  further  in  time,  until 
steady-state  conditions  were  reached,  the  surface  tenperature  will  approach  the 
steady  state  value.  This  incrH-jsed  surface  temperature  can  also  be  seen  from  Fig. 

12.  ~ 

Figure  13  shows  the  effect  of  AP-particle  size  on  radial  surface  temperature 
distribution  in  the  solid.  As  the  AP-particle  size  is  decreased,  the  radial 
tenperature  distribution  becomes  more  uniform.  This  is  reasonable  because  for  small 
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Fig.  12  Radial  Variation  of  Surface  Temperature 
Showing  the  Appearance  of  Temperature 
Peaks  near  Oxidizer -Fuel  Interface 
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Fig.  13  Effect  of  AP  Particle  Size  on  Radial  Variation  of 
Surface  Temperature  (Average  dP/dt  =120  GPa/s) 


particle  size  lateral  heat  conduction  is  more  pronounced.  In  the  limiting  case, 
as  the  AP -particle  size  becomes  vanishingly  small,  the  temperature  distribution 
will  be  uniform.  Ignition  criteria  based  upon  AP  burning  rate  (IC  2)  or  zero  gas- 
phase  velocity  (IC  3)  will  result  in  lower  ignition  delay  time  for  smaller  AP 
particle  size  because  for  smaller  particle  size  surface  temperature  on  AP  is 
higher.  On  the  other  hand,  since  the  surface  temperature  on  PBAA  is  lower  for  smaller 
AP  particle  size,  the  adjacent  gas-phase  temperature  is  lower  than  that  for  larger 
particle.  Since  pyrolyzed  oxidizer  species  can  diffuse  readily  to  the  higher 
temperature  region  over  PBAA  surface,  the  maximum  rate  of  gas-phase  heat  release 
occurs  there.  As  a  result,  for  smaller  AP  particles,  the  ignition  criterion  based 
upon  maximum  gas-phase  heat  release  (IC  1)  will  show  longer  ignition  delay  time, 
which  is  opposite  to  that  predicted  by  the  other  two  criteria  (see  Fig.  14).  How¬ 
ever,  the  predicted  ignition  delay  times  appear  to  be  datively  insensitive  to 
the  oxidizer  particle  size.  Since  predicted  ignition  delay  times  are  more  strongly 
dependent  upon  the  properties  of  fuel  and  oxidizer  as  well  as  the  chemical  kinetic 
constants,  no  conclusions  can  be  drawn  at  this  time  as  to  the  effect  of  oxidizer 
particle  size  on  ignition  delay. 


Criterion  1 
Criterion  2 
Criterion  3 


AP  Particle  Diameter,  d^p,  yin 

Fig.  14  Predicted  Effect  of  AP  Particle  Size  on  Ignition  Delays 

Using  Different  Ignition  Criteria  Based  Upon  the  Attainment 

of:  (1)  q,M  ...  . ;  (2)  r.  ...  and  (3)  zero 

M  maxjcritical  v  b, critical  v  ' 

axial  velocity  m  the  gas  phase  at,  a  specified  location. 
(Average  dP/dt  =  120  GPa/s) 
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III.  INVESTIGATION  OF  PROPELLANT  CRACK  PROPAGATION 
3 . 1  Experimental  Observations  of  Crack  Propagation  in  Burning  Solid  Propellants 

At  present,  there  is  a  lack  of  experimental  and  theoretical  knowledge  in  the 
area  of  solid  propellant  grain  fracture  during  motor  firing.  This  technological 
deficiency  becomes  more  acute  when  narrowing  the  field  to  grain  fracture  under 
rapid  pressurization  and  burning  conditions.  This  area  is  of  special  interest  to 
investigators  of  rocket  motor  grain  integrity,  where  high-speed  crack  growth  and 
propellant  fragmentation  may  provide  enough  additional  burning  surface  area  to 
enhance  the  possibility  of  deflagration-to-detonation  transition  (DDT)  in  solid 
propellant  rocket  motors. 

In  the  past,  numerous  studies,  based  purely  on  solid  mechanics  considerations, 
have  been  conducted  to  investigate  crack  propagation.  Despite  the  fact  that  there 
is  no  experimental  evidence  to  support  the  validity  of  applying  results  of  inert 
propellant  crack  propagation  studies  to  live  propellants  under  burning  conditions, 
virtually  all  previous  studies  have  employed  inert  propellants  in  non-burning 
environments.  It  is  therefore  apparent  that  an  investigation  of  crack  propagation 
under  burning  conditions  would  provide  useful  data  to  interpret  the  mechanism  of 
crack  propagation  in  a  combusting  environment . 

An  investigation  was  undertaken  to  help  bridge  this  technological  gap.  Specif! 
objectives  of  this  study  were: 

1.  To  develop  an  experimental  technique  to  study  crack  propagation  in  a 
composite  solid  propellant  under  burning  conditions. 

2.  To  measure  crack  propagation  velocities  in  solid  propellant  samples. 

3.  To  observe  any  abnormalities  of  the  propagating  crack. 

4.  To  study  the  effects  of  chamber  pressure  and  pressurization  rate  on 
crack  propagation  velocities. 

5.  To  compare  burning  crack  propagation  results  with  experimental  data  and 
theoretical  results  of  inert  propellant  crack  propagation  studies. 
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3.1.1  Experimental  Setup 

Crack  propagation  tests  were  conducted  with  a  windowed  chamber;  the  schematic 
diagram  of  the  test  rig  is  shown  in  Fig.  15.  A  propellant  sample,  with  prefab¬ 
ricated  crack,  was  internally  pressurized  by  a  small  driver  motor  which  produced 
high -temperature  and  high-pressure  gases.  The  propellant  sample  configuration 
allowed  for  significant  displacement  of  the  crack  walls,  resulting  in  high  stresses 
and  strains  at  the  crack  tip  which  in  turn  led  to  crack  propagation. 

Figure  16  presents  a  schematic  diagram  of  the  driving  motor  used  for  hot  gas 
generation.  An  electric  primer  (FA  874)  was  used  to  ignite  the  booster  propellant 
and  the  igniter  charge.  The  product  gases  flowed  through  a  multiperforated  nozzle 
into  the  main  chamber,  causing  ignition,  flame  spreading,  and  mechanical  deformation 
of  exposed  internal  crack  surfaces.  The  continued  pressurization  of  the  crack 
cavity  is  dominated  by  the  output  mast  and  energy  fluxes  of  the  driving  motor. 

Using  this  apparatus,  it  is  possible  to  obtain  pressures  up  to  50  MPa  (approxi¬ 
mately  7000  psi),  and  pressurization  rates  in  the  order  of  100  GPa/s  (approximately 
1G6  atm/s).  For  tests  conducted  in  this  study,  pressurization  rates  ranged  between 
one  and  ten  GPa/s. 

Figure  17  is  a  photograph  of  a  typical  propellant  sample,  mounted  in  a  chamber. 
The  sample  was  machined  to  a  desired  geometry  and  a  uniform  thickness.  A  3  mm 
deep  slit  was  cut  with  a  sharp  razor  blade  to  help  initiate  crack  growth.  The  pro¬ 
pellant  used  was  73/27  AP/HTPB  with  200  ym  AP.  Once  the  sample  was  installed  in 
the  chamber,  a  plexiglass  window  assembly  was  bolted  on,  slightly  compressing  the 
sample  between  the  rear  wall  of  the  chamber  and  the  window. 

A  block  diagram  of:  the  data  acquisition  system  used  in  the  study  is  shown  in 
Fig.  18.  Pressure  measurements  were  made  at  two  locations  in  the  main  chamber  (one 
at  the  crack  entrance  and  the  other  in  the  crack-tip  region)  by  piezoelectric  trans¬ 
ducers.  Signals  from  the  transducers  were  amplified  through  a  charge  amplifier 


Fig.  15  Schematic  Diagram  of  Crack  Propagation  Chamber 


Fig,  16  Schematic  Diagram  of  Driving  Motor 


Fig.  17  Photograph  of  Propellant  Sample  Mounted 
in  Crack  Propagation  Test  Chamber 


Block  Diagram  of  Data  Acquisition  System 
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before  recording  on  the  transient  waveform  recorder  end  the  high-speed  magnetic 
tape  recorder.  A  visual  record  of  the  event  was  made  with  a  16  mm  Hycam  movie 
camera,  using  a  35,000  pps  filming  rate. 

3.1.2  Experimental  Results 

Several  variables  were  considered  important  to  the  outcome  of  the  experiments: 
the  sample's  material  properties,  geometry,  and  initial  tenperature;  the  exit  nozzle 
size;  and  the  pressurization  rate  controlled  by  the  efflux  from  the  driving  motor. 

Test  results  presented  in  this  report  include  only  ef  ects  of  pressurization  rate; 
other  variables  were  held  constant. 

A  typical  pressure  vs.  time  trace  is  presented  in  Fig.  19.  From  the  plot, 
it  can  be  seen  that  during  the  time  interval  of  measured  crack  growth  (approximately 
1  msec),  the  pressurization  rate  was  approximately  constant.  Nearly  constant 
pressurization  rates  were  observed  for  all  test  firings,  and  this  fact  provided  a 
convenient  test  parameter  which  was  easily  v*  u  from  test  to  test  by  adjusting 
the  size  of  the  igniter  charge  loaded  in  the  driving  motor. 

The  bulk  of  data  obtained  in  each  test  was  recorded  by  high-speed  photography. 

A  typical  film  record  is  shown  in  Fig.  20.  Hot  gases  from  the  driving  motor  illum¬ 
inated  the  crack  cavity  and  outlined  the  crack  walls.  In  this  figure,  each  vertical 
column  of  light  is  the  outlined  crack  at  a  different  time.  The  event  proceeds  from 
left  to  right  with  a  time  interval  of  29  ys  between  consecutive  pictures. 

After  an  initial  period  of  irregular  crack-tip  displacement  (first  three  pictures) , 
the  evert  shown  in  pictures  4-9  illustrates  a  typical  mode  I  crack  propagation  with 
the  crack-tip  displacement  in  the  vertical  direction  as  expected.  However,  in  picture 
10,  a  dramatic  change  takes  place  in  the  propagation  mechanism.  A  fan  region  of 
light  emerges  just  above  the  crack  tip,  and  is  accompanied  by  a  change  in  the  geometry 
of  the  crack.  The  crack  walls,  up  to  this  point,  have  displaced  outward  towards  the 
chamber  walls.  However,  near  the  crack  tip,  a  small  triange  is  formed,  indicating 
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FIG.  19  Pressure-Time  Trace  for  Crack  Propagation  Test 
(Test  No,  DNCP1-7;  DP/Dt  =  4xl(r  MPa/s  during 

MEASURED  CRACK  PROPAGATION) 
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Fig.  20  Photograph  of  Crack  Propagation  in  a  Burning  Solid  Propellant  with 
High  Pressurization  Rate,  dP/dt  =  1.55  x  10*  at«/s  (35,000  pps) 
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the  existence  of  only  one  crack.  After  the  fan  region  appeared,  the  crack  walls 
become  almost  parallel,  and  the  crack  tip  itself  lost  its  definite  triangular 
shape  and  becaae  a  jagged  horizontal  line.  One  explanation  for  the  fan  region 
and  the  corresponding  change  in  crack  geometry  is  the  creation  of  aicropores  by 
the  material  fracture  near  the  crack  tip.  fhis  leads  to  branching  of  the  crack, 
which  allows  crack  walls  to  displace  outward  the  square-shaped  contour  described 
above . 

The  crack-tip  position  for  each  test  was  plotted  as  a  function  of  tine.  An 
example  is  presented  in  Fig,  21.  If  the  tests  were  conducted  in  a  non-burning 
environment,  one  could  easily  associate  the  crack-tip  displacement  with  crack 
propagation  velocity.  However,  in  the  case  of  a  burning  environment,  mechanical 
displacement,  material  regression  due  to  burning,  and  the  possibility  of  flame 
spreading  between  the  sample  and  window  must  all  be  considered.  Burning  rate 
calculations  reveal  that  the  effect  of  deflagration  would  only  account  for  .1% 
of  the  displacement  experienced,  and  therefore  may  be  considerably  negligible. 

The  sample's  vertical  displacement  at  the  crack  tip  due  to  pressurization  was 
calculated  using  a  finite  element  program.  From  this  analysis,  it  was  determined 
that  the  displacement  could  also  be  considered  negligible  when  compared  to  the 
total  displacement.  The  vagueness  due  to  flame  spreading  above  the  crack  tip  can 
be  eliminated  by  the  interpretation  of  the  film.  While  it  is  not  obvious  from  the 
black-and-white  reproductions  of  the  photographs  included  in  this  paper,  the  actual 
geometry  of  the  crack  cavity  can  be  identified  in  most  cases  by  using  a  motion 
analyzer  for  close  inspection  of  the  films. 

Once  the  displacement  plot  shown  in  Fig.  21  was  established  as  the  approximate 
crack  propagation  plot,  crack  velocities  were  obtained  by  a  least-square  fit.  All 
tests  revealed  nearly  constant  crack  velocities  after  a  short  initial  period  of 
irregular  crack  growth.  These  velocities  ranged  from  20  to  70  m/s.  A  major  objective 
of  this  study  was  to  determine  the  effect  of  pressurization  rate  on  crack  velocities. 
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FIG,  21  Measured  Crack  Growth  in  a  Burning  Solid  Propellant 


Since  a  constant  velocity  and  pressurization  rate  was  observed  for  each  test, 
each  experiment  provided  a  data  point  for  the  plot  presented  in  Fig.  22.  This  plot 
demonstrates  a  definite  trend,  i.e.,  higher  pressurization  rates  result  in  higher 
crack  velocities.  The  straight  line  on  the  plot  represents  the  experimental 
correlation  given  below: 

vc(m/s)  =  .169  [—■  (MPa/s)]0,608  (39) 

It  must  be  noted  that  at  present  this  correlation  has  been  developed  for  a  specific 

propellant  with  a  fixed  initial  geometry  and  temperature. 

Another  aspect  of  increasing  pressurization  rates  is  illustrated  when  comparing 

Figs.  20  and  23.  Pressurization  rates  for  the  tests  shown  in  Figs.  20  and  23  were 
5  4 

1.5  x  10  and  3.0  x  10  atm/s,  respectively.  At  the  lower  pressurization  rate,  the 
crack-tip  region  maintained  its  original  triangular  shape  with  only  small  deforma¬ 
tions,  while  the  test  conducted  at  the  higher  pressurization  rate  revealed  a  square 
crack  shape  with  a  ragged  tip  region,  indicating  crack  branching.  The  geometric 
transformation  was  accompanied  by  an  increase  in  combustion  intensity  (see  Fig.  20). 
The  pictures  on  the  right  side  of  Fig.  23  present  an  example  of  flame  spreading 
across  the  face  of  the  sample.  In  this  case,  it  is  possible  >  distinguish  between 
the  flame  spreading  and  the  triangular  crack  contour.  The  pictures  in  Fig.  20  on 
the  other  hand,  show  limited  flame  spreading  and  the  clearly  outlined  square  shaped 
crack . 

3.1.o  Discussion  of  Results 

A  major  obstacle  which  arises  when  comparing  results  of  this  study  with 
contemporary  theories  of  crack  propagation  in  viscoelastic  solids  is  the  choice 

of  correlation  parameters.  The  quasistatic  crack  propagation  model  proposed  by 

20  21  2~>  20  23 

Knauss  ’  and  Schapery,  "  and  experimentally  verified  by  Knauss,  Francis,  and 
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Langlois  and  Gonard  relates  crack  velocities  in  terms  of  stress  intensity,  Calcu¬ 
lating  the  stress  intensity  for  che  geometry  and  loading  conditions  used  in  this 
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Measured  Effect  of  Pressurization  Rate  on  Crack 
Propagation  Velocity  in  a  Burning  Solid 
Propellant  Crack 


Fig.  23  Photograph  of  Crack  Propagation  in  a  Burning  Solid  Propellant  with 
Low  Pressurization  Rate,  dP/dt  =  3  x  104  atm/s  (35,000  pps) 


study  is  complicated  and  uncertain.  The  geometry  of  the  sample  does  not  lend 

itself  to  classical  elastic  or  viscoelastic  solutions.  In  addition,  the  boundary 

conditions  change  as  the  sample’s  legs  are  compressed  against  the  chamber  walls. 

The  effect  of  the  hot  combustion  gases  at  the  crack  tip  on  the  propellant's  material 

properties  is  also  unknown.  However,  it  is  possible  to  make  a  few  realistic 

assumptions  about  the  stress  intensity  at  the  crack  tip.  The  sample  exhibits 

failure  t:  the  tip  almost  immediately  during  the  event  at  very  low  pressures  (500 

K  Pa  or  less).  As  the  load  increases,  more  energy  is  made  available  to  deform  the 

crack  walls,  therefore  increasing  the  stress  intensity  at  the  tip.  During  the 

pressurization  and  crack  propagation  interval,  the  stress  intensity  increases 

with  pressure;  however,  the  crack  propagation  velocity  quickly  reaches  a  terminal 

velocity.  This  implies  that  for  the  range  of  pressurization  rates  considered 

and  the  crack  geometry  tested,  the  crack  velocity  is  independent  of  stress  inten- 

sity.  This  phenomenon  has  been  reported  by  Kim  and  Knauss,  Swanson,  and  Gent 
27 

and  Martony  in  their  studies  of  dynamic  crack  growth.  Therefore,  it  is  believed 

that  the  present  burning  crack  propagation  study  is  in  the  dynamic  regime  and 

cannot  be  described  by  the  quasistatic  model. 

It  is  interesting  to  note  that  the  measured  burning  crack  propagation  terminal 

velocity  is  between  20  to  70  ra/s,  which  is  of  the  same  order  as  those  reported  by 
28  27 

Swanson  and  Gent  in  inert  propellants  and  rubber,  respectively.  Swanson 
reported  maximum  velocities  of  approximately  40  m/s,  while  Gent  ob.*erved  crack 
velocities  of  80  m/s  in  specimens  with  no  pre-imposed  strains  in  the  crack  direc¬ 
tion.  Both  Swanson  and  Gent  proposed  upper  bounds  to  dynamic  crack  propagation 
rates  based  upon  characteristic  wave  speeds  in  the  material.  Swanson  selected 
the  glassy  Ra/leigh  wave  velocity  as  an  upper  bound  to  crack  propagation  velocity. 

He  found  that  for  inert  propellants  the  maximum  velocities  measured  were  one  order 
of  magnitude  below  th<->  Raleigh  wave  speed.  Gent  proposed  that  the  crack  propagation 
velocity  in  rubber  was  proportional  to  the  velocity  of  sound  in  the  material. 


More  specifically,  the  crack  propagation  velocity  V£  is  related  to  the  speed  of 
sound  Vs  by: 

Vc  -  .3  Vs  (40) 

Gent  also  noted  that  the  speed  of  sound  is  strongly  dependent  on  the  state  of 
strain  in  rubber. 

Dynamic  material  properties  for  the  propellant  used  in  this  study  are  not 
fully  defined;  however,  typical  values  for  similar  propellants  were  assumed,  and 
sample  calculations  were  performed  in  order  that  comparisons  could  be  drawn.  One 
can  approximate  the  speed  of  sound  in  a  viscoelastic  material  using  the  formula 

> 

given  below: 

Vg  =  v€7p  (41) 

where  E  is  the  dynamic  tensile  modulus  and  p  is  the  density  of  the  material.  Using 
Eq. (41) ,  the  calculated  Vg  is  20S  m/s  for  a  dynamic  tensile  modulus  of  63  MPa 
and  a  propellant  density  of  1.49  gm/cc.  E  was  determined  from  an  experimentally 
derived  plot  of  E  as  a  function  of  loading  time,  and  a  maximum  loading  rate  of 
24,000  MPa/s  was  used.  Using  Eq.  (40),  is  found  to  be  62  m/s.  The  actual 

crack  propagation  velocity  measured  for  the  maximum  pressurization  case  was  68 
m/s,  which  is  close  to  the  calculated  value  obtained  from  Eqs.  (40)  and  (41). 

Even  though  there  is  some  agreement,  the  method  of  calculating  the  speed  of  sound 
is  crude,  since  V  depends  on  the  state  of  strain  in  the  sample,  as  pointed  out  by 
Gent,  and  may  also  be  affected  by  the  geometry  of  the  sample  and  the  hot  gas  envir¬ 
onment  . 

At  present,  the  strong  dependence  of  crack  propagation  velocity  on  pressurization 
rate  (as  illustrated  by  the  correlation  given  in  Eq.  39),  is  not  fully  understood. 

One  explanation  for  this  dependence  is  the  transition  of  the  material  towards  a 
more  brittle  state  as  loading  rate  i3  increased.  Using  Eqs.  (40)  and  (41)  and  an 
experimentally  derived  expression  relating  the  dynamic  tensile  modulus  to  pressuri¬ 
zation  rate,  one  can  examine  the  approximate  effect  of  the  material  transition  on 


c.  Two  provide  for  a  longer  period  of  crack  growth  and  therefore 
increase  the  quantity  of  data  recorded  for  each  test. 

3.2.2  Increased  Chamber  Depth 

a.  To  house  the  igniter  assembly  in  the  main  body  of  the  crack  propagation 
chamber.  By  so  doing  the  hot  gas  passage  is  located  away  from  the 
interface  between  main  body  and  window  retainer  piece;  and  therefore 
any  possible  leaks  can  be  prevented. 

b.  To  allow  for  thicker  plexiglass  windows  for  strengthening  the  test 
chamber  and  also  to  reduce  the  possibility  of  flame  penetration  between 
the  propellant  sample  and  sacrificial  window. 

3.2.3  New  Instrumentation  and  Lighting 

a.  Two  different  locations  of  feedthroughs  for  mounting  break  wires 
or  thermocouples  to  detect  the  arrival  of  hot  combustion  product 
gases  and  instantaneous  crack  tip  front. 

b.  Additional  pressure  transducer  ports  above  the  initial  crack  tip. 

c.  External  lighting  by  using  a  light  source  with  fiber  optics  to 
illuminate  the  sample  surfaces  not  in  contact  with  combustion  gases. 
This  provides  informat ion  on  the  outer  boundaries  of  the  propellant 
sample. 

d.  Gas  sampling  ports  are  provided  for  collecting  and  analyzing  combustion 
product  gases. 

3.2.4  Strength  Considerations 

A  detailed  stress  analysis  was  done  for  a  chamber  pressure  of  5000  psi  to 
determine  the  chamber  wall  thickness,  plexiglass  window  thickness  and  the 
required  number  and  size  of  bolts  for  retaining  the  window  assembly,  igniter 
system  and  the  exit  port. 


so 


A  thick  (3-1/2")  window  retainer  plate  with  two  separate  viewing  areas, 
instead  of  a  single  long  opening,  was  designed  in  order  to  prevent  window 
buckling  and  retainer  plate  warpage.  The  front  views  of  the  window  retainer 
plate  is  shown  in  Fig.  25. 

3.2.5  Test  Chamber  Fabrication 

304  stainless  steel  material  was  chosen  for  the  test  chamber  fabrication 
in  order  to  resist  corrosive  effects  of  combustion  product  gases.  The  test 
chamber  components  are  being  constructed  at  Houtz  Instrument  and  Machine 
according  to  detailed  mechanical  drawings  and  the  fabrication  is  expected  to 
be  finished  by  November  30,  1982. 

3.3  Finite  Element  Analysis  of  Propellant  Sample  Used  in  Crack 
Propagation  Studies 

A  stress  analysis  of  the  propellant  sample  was  performed  in  order  to  gain 

some  insight  into  the  mechanical  behavior  of  the  sample  material  under  constant 

/ 

internal  pressurization  of  the  crack  cavity.  Several  factors  make  this  analysis 
complicated  and  therefore  certain  assvmiptions  were  made  to  simplify  the  problem  and 

reduce  computation  time.  These  assumptions  reduce  the  accuracy  of  the  results; 

/ 

however  the  analysis  provides  a  qualitative  picture  of  the  strain  and  stress 

/ 

distributions  in  the  sample  before  crack  propagation. 

The  finite  element  grid  for  the  model  was  generated  on  an  interactive  mesh 
generation  program  (STAB) .  The  sample  was  considered  to  be  symmetric  about  the 
crack  axis  and  a  finite  crack  radius  of  1.78  mm  (.07  in.)  was  used.  The  boundary 
conditions  on  the  crack  sample  for  the  finite  element  program  are  shown  in  Fig.  26. 

The  model  was  solved  using  a  non-linear  finite  element  code  (ABAQUS) .  The 
assumptions  and  options  used  are  listed  below: 

1.  An  elastic  analysis  was  used;  however,  the  sample's  viscoelastic  behavior 
was  approximated  by  using  the  relaxation  modulus,  corresponding  to  the 


E  *  8000  psi 
v  ”*  500 
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test  pressurization  rate  and  initial  tenperature,  in  place  of 
the  elastic  tensile  modulus.  This  is  a  comoti  procedure  used  in 
solid  propellant  grain  deformation  analysis. 

2.  The  program  allowed  for  incompressible  materials  end  a  realistic 
poisson  ratio  of  .5  was  used  for  the  propellant  sample. 

3.  A  non-linear  large  deformation  analysis  was  performed.  This  option 
is  needed  since  the  observed  mechanical  deformation  is  quite  signifi¬ 
cant.. 

4.  Non-linear  boundary  conditions  were  modeled  using  gap  elements  between 
the  sample  and  the  chamber  wall.  This  allows  for  unrestricted  displace¬ 
ment  of  the  sample  lags  until  the  outer  surface  is  in  contact  with  the 
chamber  wall. 

5.  This  preliminary  analysis  did  not  allow  for  the  propagation  of  the 
crack  tip  since  ABAQUS  does  not  have  this  capability.  In  spite  of  this 
fact  the  solution  obtained  are  extremely  useful  for  the  understanding 
of  the  mechanical  behavior  of  the  sample  prior  to  fracture. 

The  finite-element  stress  analysis  results  are  presented  in  Figs.  27  and  28. 
Figure  27  is  a  displacement  plot  where  the  dotted  lines  represent  the  initial 
location  of  the  sample  at  initial  pressure  of  0.1  MPa  and  the  solid  lines  repre¬ 
sent  the  displaced  elements  at  5.17  MPa.  It  can  be  seen  clearly  from  this  figure 
that  the  mechanical  displacement  of  the  crack  contour  is  significant;  however,  the 
initial  triangular  shape  of  the  crack  cavity  is  maintained.  Figure  28  shows  the 
contours  of  maximum  principal  deviatoric  stress  near  the  crack  tip  region.  An 
interesting  fact  to  note  from  this  plot  is  that  the  maximum  stress  is  located  near 
the  shoulder  of  the  crack  tip  instead  of  at  the  apex  of  the  crack.  This  is  due 
to  large  deformations  perpendicular  to  the  crack  axis.  The  result  supports  the 
observed  transformation  of  the  crack  contour  from  initial  triangular  into  a  square 
shaped  contour  as  discussed  in  Section  3.1- 
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Fig.  27  Calculated  Displacement  of  the  Propellant  Crack 
Sample  from  Finite  Element  Program 
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Fig.  28  Calculated  Contours  of  Maxiaura  Principal  Deviatoric 
Stress  Near  the  Crack  Tip  Region 


IV.  SUMMARY  AND  CONCLUSIONS 


Some  important  observations  and  conclusions  from  the  crack  tip  ignition  study 
are  summarized  as  follows: 

1.  A  bright  luminous  zone  behind  the  reflected  compression  wave  was  observed 
at  high  pressurization  rates.  Diagnostic  experiments  reveal  that  the 
luminous  zone  is  caused  by  combustion  of  unreacted  species  from  the 
igniter  system. 

2.  Experimental  results  indicate  that  as  the  pressurization  rate  is  increased, 
heat  flux  to  the  propellant  surface  increases  and  hence  ignition-delay  time 
decreases.  No  distinguishable  effect  of  the  crack  gap  width  on  the  ignition 
process  was  evident  for  the  range  considered.  Limited  results  obtained  by 
using  an  aluminized  propellant  as  the  igniter  show  that  the  ignition-delay 
time  is  somewhat  lower  for  higher  flame-temperature  igniter  gases. 

3.  Calculated  ignition  delay  tj.mes  based  on  several  criteria  (namely,  attain¬ 
ment  of  a  critical  volumetric  gas-phase  heat  release  rate,  a  critical 
pyrolysis  rate,  and  zero  axial  velocity  in  the  gas-phase  near  the  propel¬ 
lant  surface)  are  in  agreement  with  the  experimentally  observed  trend  that 
Ignition  delay  time  decreases  with  increasing  pressurization  rate.  Pre¬ 
dicted  values  of  ignition  delay  for  lower  pressurization  rates  showed 
greater  dependence  on  the  choice  of  ignition  criterion. 

4.  The  solutions  revealed  that  the  ignition  process  consists  of  the  following 
sequence  of  events:  (a)  As  pressurization  begins,  the  gas-phase  tempera¬ 
ture  starts  to  rise,  the  unburned  species  near  the  propellant  surface 

react,  and  heat  is  transferred  to  the  propellant  surface,  (b)  Following  a  per 
iod  of  inert  heating  and  continued  pressurization,  the  surface  temperature 
rises,  and  oxidizer  and  fuel  species  pyrolyze.  (c)  When  the  concentration 
of  o.ddizer  and  fuel  species  and  the  local  gas-phase  temperature  are 
sufficiently  high,  intense  gas-phase  reactions  begin. 
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5.  The  augmentation  in  heff.t  feedback  to  the  propellant  at  higher  pressurization 
is  a  result  of  a  combination  of  the  following  mechanisms:  heating  due 

to  compression  wave  reflection  at  the  closed  end;  heat  release  due  to 
burning  of  unreacted  igniter  species  near  the  tip,  behind  the  reflected 
compression  wave;  and  enhanced  heat  transfer  due  to  recirculating  hot  gas 
near  the  tip. 

6.  Solutions  revealed  significant  temperature  non -uniformity  implying  two- 
dimensional  effects  which  must  be  included  in  modeling  composite  propellant 
ignition. 

7.  More  accurate  chemical  kinetic  constants  pertaining  to  ignition  are  needed 

to  describe  the  pyrolysis  and  gas-phase  reaction:,  for  more  precise  prediction 
of  ignition  delay  times. 

Important  conclusions  that  can  be  drawn  from  the  crack  propagation  investigation 
are  described  below.  Even  though  the  study  presented  in  this  report  is  only  a  first 
step  in  the  investigation  of  crack  propagation  in  a  burning  solid  propellant,  sev¬ 
eral  interesting  facts  about  the  phenomenon  have  been  observed  and  are  summarized 
below: 

lt  Crack  propagation  velocities  of  20  to  70  m/s  were  measured  in  a  burning 
composite  solid  propellant.  These  velocities  are  in  the  same  range  as 
those  reported  in  non -burning  dynamic  testing  of  inert  propellants  and 
rubbers . 

2.  It  was  found  that  cracks  propagated  at  a  constant  velocity  when  subjected 
to  a  constant  rate  of  internal  pressurization. 

3.  It  was  observed  that  there  is  a  characteristic  difference  in  crack  propaga¬ 
tion  processes  between  low  and  high  rates  of  pressurization.  The  crack 
lip  geometry  maintained  its  original  contour  for  low  dp/dt .  However,  a 
fan  region  above  the  crack  tip  (accompanied  by  a  squaring  off  of  the  tip 
region)  was  observed  for  high  pressurization  rates  in  the  order  of  20,000 
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MPa/s.  The  uneven  multiple  front  of  the  crack  tip  and  the  appearance 
of  the  fan  region  suggests  the  occurrence  of  crack  branching  and  micro- 
structure  damage. 

4.  A  tentative  experimental  correlation  showed  a  stronger  dependence  of 

crack  propagation  velocity  on  pressurization  rate  than  predictions  based 
solely  on  the  elastic-glassy  transition  behavior  of  viscoelastic  material 
during  dynamic  loading.  The  reason  for  this  discrepancy  is  unknown,  and 
further  research  in  the  propagation  of  burning  solid  propellant  cracks  is 


needed. 
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APPENDIX  A 


Details  of  the  Numerical  Scheme 


The  coordinate  transformation  described  by  Eq.  (37)  and  Eq.  (38)  is  used  to 
map  the  (r,  z)  coordinates  in  the  physical  space  to  (r,  s)  coordinate  in  the  trans¬ 
formed  space.  The  transformed  governing  equations  for  the  solid  and  gas  phases 
take  the  following  form. 

(A)  Solid  Phase  Equations: 


Oxidizer: 

pox.5Co*.5  W  ‘  VstA>!>2  0  *  £> 


(A-l) 


Fuel: 

n  r  £L 

PF,SLF,S  3t 


(1*s)2$aa' 


kP,S[7  37*01^F,s 


(A-2) 


(B)  Gas  Phase  Equations: 
Continuity: 


V  a 

H4*  V*-s) 


3(PgVz) 


(A- 3) 


Energy: 

cppg  w *  'Ypvz  -  H 


*  V*(1-)2  0 *  **  E  *  ^r>  ‘  <g’ 


(A-4) 


Species 


8Y.  9Y.  3Y. 

pg  5T-  *  p8Vg(1-s)  rr  *  f  If  <rD»g  sr*  *  “i 


CA-5) 
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where  i  «  1,  2,  3,  and  4  represent  oxidizer,  NHj,  HCIO^,  and  fuel,  respectively. 
A  second  transformation  is  used  to  map  the  (r*  s)  coordinates  into  (r,  s'). 

This  transformation  is  described  by  the  following  equation, 
s'  ■  -s  for  -1  <  s  <  0 

(A-6) 

s'  «  s  for  0  £  s  £  1 

This  transformation  affects  only  the  solid  phase  equations  and  these  equations  take 


the  following  form. 


Oxidizer: 


c  LL 


OX, S  OX, 5  3t  OX,Sl  S 


Kox,s[Ast1-S,)2ri’Asf1-S,)lfl 


+  K-  f—  |X  +  Li  +  qti 
Ox,slr  3r  9r 2J  4ox,s 


(A-7) 


Fuel: 


PF.sS.sfr^F.c^^’^ 


Hi  -  ivH 


+  k  rJL  LL  +  LII]  +  A.. 

ox,slr  3r  dJ,2l  qF,< 


(A-8) 


The  gas  phase  equations  are  obtained  by  replacing  s  by  s'  and  are  mathematically 
identical  to  A-3,  A-4  and  A-5. 

F.D.  Approximations  for  Radial  Derivatives  Using  Allen's  Method 

A  three  point,  variable  mesh,  Allen's  Method,  is  used  to  approximate  derivatives 
in  the  radial  direction.  Allen's  Method  is  chosen  in  preference  to  Central  Difference 
Method  because  it  offers  a  better  approximation  of  the  derivative  near  the  axis  of 
the  computational  domain.  Figure  A-l  shows  a  z  =  constant  plane  near  the  axis 
showing  nodes  with  different  radii. 

The  following  paragraphs  describe  the  derivation  of  the  difference  formula  using 


Allen's  Method.  The  radial  derivative  term  appears  in  the  following  form  in  our 
governing  equations. 
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3D 


1  ■ 


s  ■ 


e 


•— -y  +  ■—  (where  n  could  be  temperature  or  species  mass  fraction) 
8r 

We  define  the  following  quantities: 

n  =  I|l 

9r2  r5F 

ft  is  assumed  to  be  locally  constant 

Hum  JF  *  F  r  "  8 


The  solution  for  Eq.  (A-ll)  is 


••• n  ■  l 


r  +  ~ 
r 


Tdr  *  j  r2  +  A  Jtnr  +  B 


(A-9) 

(A- 10) 
(A-ll) 

(A- 12) 

(A- 13) 


To  evaluate  the  constants  ft.  A,  B  in  Eq.(A-13),  we  write  rj  at  three  consecutive  points 

j-1.  j+l- 


Vi 

ft  2 

"  4  Tj-1 

+  A  tnr^  1  ♦  B 

(A-14) 

t 

ni 

*  j  r2  + 
4  ] 

A  Jt  nr.  +  B 

J 

(A-15) 

Vl 

=  r2 

4  j+1 

+  A  *-nrj+1  ♦  B 

(A-16) 

We  solve  Eqs.  (A-14) ,  (A-15),  and  (A-16)  to  get  the  following  value  for  ft. 


Treatment  of  Nodes  on  the  Axis  and  One  Node  Away  from  the  Axis 
(a)  Approximation  of  Radial  Derivative  Term  on  the  Axis : 

Assume  that  the  scalar  variable n  is  symmetric  about  the  axis;  thus, 


4  tnCrj+i/rj5  4£n(*/r.  ) 

0 - <Vi-v  *  — H4-  (viV 

(A-17) 

where  Q  *  An  (r^^)  (r^  j-r2)  +  in  (r^/ri  _l)  (r|+1-r2) 

(A- 18) 

.-i,,-  u*.  a  * 


Use  the  Taylor’s  series  to  expand  the  tern  ^  about  r-o 


15-  'jfr*  ”  *  r  A-  >  ♦ 

9r  r*o  3r“  r»o  2  sr3 


3r  r*o 


We  neglect  the  terms  with  order  higher  than  2,  and  divide  Eq.  (A-20)  by  r 


I  In  -  'o.  ) 
r  3r~i7}r.o 


Thus  for  the  nodes  on  the  axis 


I  |n.  +  .  2  £_i 

r  3r  23r2  3? 

We  define  ■  fl’ 

3r 


From  Eqs.  (A-22)  and  (A-23) 


IS*n'r 

n  *  y-r2  +  B’ 

Qt  2 

n  *  y-  *  ♦  B' 
at  r  *  o  »  r. 


nx  -  B' 

n  2  ~  y  1*2  *  r)  j 

2(n,-n,) 

Thus,  fi'  *  — — —— 
r2 


From  Eos.  (A-22),  (A-23),  and  (A-25) ,  we  have 

l  3n  ,  32n  _  ,  32n  _  4(rW 
r  Sr  2  - * - 


3r 


for  nodes  on  the  axis. 


(b)  Approximation  of  the  Radial  Derivative  Term  One  Node  Away  from  the  Axis 

Consider  the  most  general  expression  for  n  given  by  Eq.  (A- 13).  As  we  approach 
the  axis,  we  have 

n  *  Him  (^-  r2  ♦  A  Jlnr  +  B) 
m> 

*  -®A  +  B 

This  result  is  unreasonable,  hence  we  conclude  that  A  ■  0  near  the  axis. 

Thus,  very  near  the  axis  we  have 

n  *  j  r2  +  B  (A-27) 


We  use  Eq.  (A- 27)  for  node  2  and  3 


n2  *  I  r2  +  B 
n3  *  f  r3  +  B 


Solving  Eqs.  (A-28)  and  (A-29)  simultaneously  we  get 

4<W 


n 


r2-r2 
3  r2 


Substituting  this  expression  in  Eq.  (A-9) ,  we  get  for  node  2 


3jn  .  1  3n  4(n3‘n2i 

Z2  r  -5F  *  -t;2 

3  r3  r2 


(A-28) 

(A-29) 


(A-30) 


(A-31) 


Solution  of  the  F.D.  Equations 

A  generalized  Crank  Nicolson  Algorithm  was  used  to  solve  the  finite  difference 
equations  obtained  by  the  procedure  described  above.  Figure  A-2  represents  the  radial, 
axial  and  time  coordinates  for  the  above  scheme.  The  algorithm  developed  had  a 
facility  to  make  the  method  either  fully  implicit  or  semi-implicit  in  time. 

Following  equations  are  obtained  for  the  derivative  terms  in  the  governing 
equations.  The  parameter  8  determines  the  implicitness  of  the  method.  A  central 
difference  scheme  is  used  to  approximate  the  axial  derivatives,  while  Allen's  method 
described  in  the  previous  section  is  used  for  the  radial  derivatives. 


Temporal  derivative 

j  k+6 

3n  &n 


Axial  derivative 


Radial  derivative: 


l  3n  „  3‘n  1  3r, 


!  k+e  9  !k+e 

7  i 


£1L  +  _  'I  a  ±  ZlLI  +  1 

3r  3r2  r  5r|  i  j  3l>2  ^ 


For  internal  nodes,  the  radial  derivative  can  be  expressed  as 
,  ^2_  4  An(r..,/r,)  ^  ,  4  Jnfr./r, 


1  3n  A  3  l 
?5r%; 


,2  1  Q  Wi,J-l  ni,j’  <1  Wi,3-* 


*  (1-9)  [ 


4  *n(,Wrj)  ,  It 


On?  .  ,  -  n?  •)  + 
i.j-i  i>r 


4  hifrj/Tj.,)  t 
Q  tni 


(\.m  ■  ni 


where  Q  is  defined  l y  Eq  (A- 18). 

For  nodes  on  the  axis  and  one  node  away  from  the  axis  (r  =  o  =  and  r 
at  r  *  o  *  rj  (node  on  the  axis) 


i  ,  „k+l  ,  k+1  ,  ^k  .  k 

r  57  *  — T  ■  e( - ^7 - ^  *  (1-9)  ( - -^—2 - £ii) 

3r  r2  r2 


r  37 


at  r  *  r2  (first  node  away  from  the  axis) 

i  a  a2  4  ~  4  n*+l  4  n*  , 

-  8  (—*4 - 7-^-)  *  (1-6)  (~i£ 


.  k 
4  n.  ^ 
_ x_» -) 


Quasilinearization  of  Source  Terms 

The  source  terms  in  energy  equations  and  species  equations  are  functions  of 
temperature  and  species  mass  fraction.  These  quantities  being  the  unknowns 
for  which  energy  and  species  conservation  equations  are  solved,  quasi lineariza¬ 
tion  of  the  source  terms  is  required. 

Let  F(T,Yj)  represent  the  general  inhomogeneous  teim.  By  Taylor  series 
expansion  we  write 


pK+8 

Fl+1 


„K+6 


fTK*6 

UJU1 


-Tf9) 


3F 

■ST 


K+8 


C* 


K+e 

j,£+i 


vK+e. 

"  Yj,*J 


3F 

577 

3 


K+e 


(A- 38) 


The  index  K+6  represents  a  time  step  in  the  Crank-Nicolson  algorithm,  where  for 
any  variable  U, 


UK+6  «  0UK+1  +  (1~8)UK 


(A- 39) 


The  index  i  represents  the  number  of  iterations  of  the  quasi linearization  loop, 

and  j  denotes  different  species. 

3F  3F 

The  derivations  -gy  and  -gy-  are  evaluated  from  the  chemical  reaction  modelling 

3 

described  earlier. 

The  procedures  described  in  this  Appendix  are  incorporated  in  the  program 


outlined  in  Fig.  3. 


Fig.  A-l  Two-Dimensional  Finite  Difference 
Cylindrical  Coordinate  System  Near 
the  Axis  Showing  a  Z* Const ant  Surface 


Fig.  A-2  Generalized  Crank -Nicolson  Coordinates 
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